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Relevance 



Mathematical modeling of many physical processes such as diffusion, vis- 
cosity of fluids and combustion involves differential equations with small 
coefficients of higher derivatives. These may be small diffusion coefficients 
for modeling the spreading of impurities, small coefficients of viscosity in 
fluid flow simulation etc. 

The difficulty with solving such problem is that if you set the small param- 
eter at higher derivatives to zero, the solution of the degenerate problem 
doesn't correctly approximate the original problem, even if the small pa- 
rameter approaches zero; the solution of the original problem exhibits the 
emergency of a boundary layer. As a result, the application of classical 
difference schemes for solving such equations produces great inaccuracies. 
Therefore, numerical solution of differential equations with small coefficients 
at higher derivatives demands special difference schemes exhibiting uniform 
convergence with respect to the small parameters involved. 

Current status of this problem 

The development of numerical methods for problems with exponential and 
power-law boundary layers introduced among others the following methods: 

i) densification of grids in boundary layers; 

ii) scheme adjustment on the boundary layer component of solutions. 

The first approach includes the works of N. S. Bakhvalov, V. D. Liseikin, 
G. I. Shishkin, V. B. Andreyev, R. Vulanovich and other authors. In 
P. S. Bakhvalov's approach, the grid nodes distribution is constructed in 
such a way that approximation inaccuracy on the grid nodes of the bound- 
ary layers stays equal while outside of the boundary layer the mesh remains 
uniform. 

It is shown that application of such a grid yields second-order accuracy with 
respect to the number of grid points. V. D. Liseikin purposes to perform 
change of variables in such a way that derivatives up to some order are 
uniformly bounded. In terms of original variables, this is equivalent to mesh 
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densification. G. I. Shishkin defined an approach for building grids that are 
uniform both within the boundary layer and outside of the boundary layer. 
In the works of G. I. Shishkin, V. B. Andreev, E. A. Savin, N. V. Kopteva 
it is shown that on such a grid, a number of difference schemes (including 
non-monotone scheme of central differences) are uniformly convergent. The 
significance of this approach for partial differential equations is strengthened 
by the fact that in the case of a parabolic boundary layer, as shown by 
G. I. Shishkin, there is no uniformly convergent adjustment schemes on a 
uniform grid. 

The second approach is presented in the works of A. M. Ilyin, G. I. Shishkin, 
K. V. Emelyanov, D. Miller, R. Kellogg and others. The main idea of this 
approach is isolating the boundary layer component of solutions and con- 
structing a difference scheme which is exact on the boundary layer function. 
The advantage of this approach is that it does not impose restrictions on the 
mesh gauge, while the drawback is that the boundary layer function must 
be expressed explicitly and the scheme must be specifically adjusted to the 
function, which is not always possible. 

In this article author investigates two nonlinear boundary value problems 
on a finite interval, resulting in exponential and power-law boundary layers. 

In the first chapter author constructs a first-order accurate finite difference 
scheme for an exponential boundary layer problem. Author then introduces 
a certain grid where on each interval the coefficients are replaced by con- 
stants. The second chapter provides a construction of a first-order accurate 
difference scheme for a problem with power-law boundary layer using a sim- 
ilar technique with additional mesh densification on the boundary layer for 
scheme convergence. 
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Chapter 1 



Nonlinear differential 
equation with exponential 
boundary layer 



The purpose of this chapter is to construct a difference scheme for a nonlinear 
equation with exponential boundary layer, uniformly convergent in the small 
parameter. Since the presence of a small parameter e at the second-order 
derivative causes an exponential growth of solutions near the boundaries 
of the interval, application of classical difference schemes for solution of 
equations with small coefficient at the highest derivative produces great 
inaccuracies. This follows from the fact that the accuracy of known schemes 
is estimated as a product of some power of the mesh size h and maximum 
absolute value of some derivative of the solution, and derivatives of the 
solutions grow unbounded in boundary layer with decreasing e. Thus the 
classical accuracy estimation of difference schemes for solving equations with 
small parameter is not acceptable. 

Let us consider the boundary value problem: 



Let's assume the function g(s) to be twice continuously differentiable for all 
s G R, e G (0, 1], and a(x) continuously differentiable. 




T £ u = -eu" + a(x) v! + g(u) = 
u(0) = A, R £ u = eu(L) + f(u(L)) = 




assuming that: 




(1.0.2) 
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1.1. ANALYSIS OF THE ORIGINAL PROBLEM 



C and Ci will everywhere designate positive constants, not dependent on e 
and the mesh gauge. 

1.1 Analysis of the original problem 

Define an auxiliary linear operator L: 

Lu = —eu" + a(x) u + c(x) u 
with boundary conditions: 

tt(0), Du = eu'(L) + d(x)u. 

Let us find out whether for the operator L the Principle of Maximum holds. 
Lemma 1.1.1. Assume 3(j)(x) > : L(p(x) > 0, Deft > 0. Then for the 
operator L the Principle of Maximum holds, i.e. for an arbitrary twice 
differentiable function ip{x) the condition: 

Lip(x) > 0, D^j > 0, ip(0) > 

implies that: ip{x) > 0. 

Proof. Assume i(j(x) to be some differentiable function. Express ip(x) as the 
product: ip(x) = v(x)cf>(x). 
Assume there is xo : ip(xo) < 0. 

Then v(xq) < 0. 

v(0) > 0, Dip = {y4>)\L) + dv(j)(L) > 
v'4> + v4>' + dvcj) > 0. 
v'cj) + [<f/ + d(f)]v > 0. 

Consequently, 

v'<j) + D(j)v > 0. 

Lip = L(v(p) = -s(v(p)" + a(v(p)' + cv(j) 
= -e(v'(j> + V(j)')' + a(v'4> + vcj)') + evej) 
= -e{v"(j) + 2v'4>' + vcj)") + a(v'cf) + vcf)') + evej) 
= -e(v"<t> + 2v'0') + v(-e<f>" + a<f> + c(f>) + av'cj) 
= -ev"<j) + 2ei;V / + av'cj) + «L0 

= -e(j)v" - v'(2e<t>' - a<t>) + vL(p > (1.1.1) 
Consider two cases: 
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1.1. ANALYSIS OF THE ORIGINAL PROBLEM 



i) Let v(L) > 0. Then there is 77, such that v(rj) < holds on the interval 
[0, L\. Then rj is a minimum point for the function v(x). At the point 
of minimum first derivative equals zero, therefore, given the conditions 
of Lemma, in (1.1.1) £(j)v" > 0, vLcj) > 0, and the term v'{2e(j)' — a<p) 
vanishes. So Lip < 0, but from condition of the Lemma Lip > 0, a 
contradiction follows. 

ii) Let v(L) < 0. Then due to the inequality 

v'cj) + D(j)v > 

we obtain 

v'(L) > 

and, consequently, there is a minimum point fj. Arguing similarly to 
the first case, we arrive at a contradiction. 

□ 



Suppose that 

d(x) > 0, c(x) > -j3, a 2 - 4/3e > 7 > 0, j3 > 0. 

We show that for L the Principle of Maximum holds. Define the function 

(f){x) as: 

Mx) = e l ix ~ L \ where 5 = —. 

a 



Dcj)(L) = 5 + d>0. 



> 



(- 


8 2 




e 


(- 


6 2 




e 






^2 


(a 2 



a- 



Thus, 

<t>>0, D<p>0, L(f)> 0, 
i.e. the Principle of Maximum holds. 

Now we obtain an estimate for stability for the problem (1.0.1). 
Lemma 1.1.2. Let p(x) and q(x) be two arbitrary functions. Then 

\\p(x) - q(x) || < c\\T £ p-T e q\\ + c\p(Q) - q(0) \ + c\R e p - R e q\ 
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1.1. ANALYSIS OF THE ORIGINAL PROBLEM 



Proof. Define z = p — q. 

Then with respect to z we obtain the boundary problem: 

Lz = —ez" + a(x) z + c(x) z = T e p — T £ q, where 

c{x) = 9{P) ~ 9{q) > -3. 
p-q 

z(0)=p(0)-g(0). 



Dz = ez'(L) + dz(L) = ez'(L) + ~ f ^ z(L) 

p{L) - q{L) 

= ez'(L) + f\0) = R £ p - R e q. 
Introduce the function 

$( x ) = Cie ^ ||T e p - T s q\\+c 2 e^ x -V \\R eP - R £ q\\+e^ \p(0) - q(0)\±z(x) 

It is easy to verify that: 

„2 „,2 \ / „,2 



a(x — L) 

Le 2i > 



f a 2 -40e \ ea{x - L ) 2£ > 7e a(,-i) 2£ > 



Le ^ > + 2(3 + c (x)] e 2 ^a > f-e^ -f 



/3e^a ( > > 0. 

ar / ar 



Consequently, 



14 > C!^e^ x a \\T £ p - T £ q\\ - \\T £ p - T £ q\\ > 0, 



where c\ 



a 

„2 



Pi' 



28x / 28 \ 2Bx 

Be~ = ( e— + d ) e~ > 0, 

a 

a.(x — L) 

De 2 E = a + d > 0. 
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1.2. ESTIMATE OF THE DERIVATIVE 



Consequently, 



Lip > (a + d)c2 \R £ p — ReQ\ — \ReP — R e q\ > 0, 

where c 2 = 
Thus, defining 

Vw = — 2 e a \\ T eP-T £ q H — e * ±z(x), 

we obtain ip(0) > 0, Dtp > 0, L?/> > 0, hence by Lemma 1.1.1 ip{x) > 0. □ 

Corollary 1.1.3. Lemma 1.1.2 implies uniqueness and boundedness of prob- 
lem (1.0.1). 

Proof of the boundedness: 

Let p(x) = u(x), q(x) = 0, then by Lemma 1.1.2, 

O? 2f3x 1 a(x-L) 2/3x 

\u(x)\ < —e— \g(0) + — -je— i— |/(0)| + e" |u(0)| . 
The uniqueness of this solution is obvious. 

1.2 Estimate of the derivative 



We obtain an estimate of the derivative problem (1.0.1). 
Lemma 1.2.1. 

i 1 1 \ i 1 a ( j: - L ) 

p {x)\ < c 



l + -e" 

e 



Proof. We express the equation (1.0.1) in the form: 

' r r L 

+ g(u(x)) exp 



^— eu' exp 






7x e 



a(t) 



dt 



Integrating from x to L: 
—eu'(L) + eu'(x) exp 

u'(x) = u'(L) exp 





hi 















a(t) 



r/f 



s(u(s))exp 



(is = 



o(t) 



dt 



ds 
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1.3. CONSTRUCTION OF A DIFFERENCE SCHEME 



It is clear that: 



g(u(s))exp 



a(t) 



dt 



ds < - I 



a(t) 



dt 



ds 



ef^-V ds 



£ 
C 

< - = Cl 
a 



From the boundary condition we conclude that 

f«L)) 



\u'(x)\ < 



c 

< - 



consequently, 



\u'(x)\ <c+ -£^- L \ 



□ 



The obtained estimate of the derivative characterizes the boundary layer. 

1.3 Construction of a difference scheme 

Let us introduce a non-uniform grid and replace the coefficients by constants 
in each interval of that grid. That allows us to write the solution explicitly. 
Matching derivatives on adjacent intervals will lead to a difference scheme. 

So, let A n = [x n -i,x n ]. 

Now we turn to a problem with piecewise constant coefficients: 

j eV" - aV' + g(V) = 0, 
{ V(0) = A, R e V = 0, 

where a = a n = a(x n -i), g(V(x)) = g n = g(V(x n -i)), at x G A n . 

On an arbitrary interval A n we have: 

J eV" - a n V + g n = 0, x E A n 
[V(x n - 1 ) = Vti, V(x n ) = V^ 



1.3. CONSTRUCTION OF A DIFFERENCE SCHEME 



where {V^ 1 } is not yet denned. 

The solution on the interval A n has the form: 

V(x) = ci + c 2 e^ {x - Xn) + —x. (1.3.1) 

Let h n be the grid size: h n = x n — x n -\. 
Define c 2 from the boundary conditions: 

' Cl+C2e -^ n + 9n Xn _ i= yH^ 

ci + c 2 h — £ n = 



Let us find c 2 : 



'-X n -1 = VI 



C2 X n , 



v n-l v n + „ n 



a,, 



C2 e 



C2 



1 • 



Then let us find V'(x): 



e a 



n 



To ensure that the solution is continuously differentiable, the derivatives of 
solutions must be matched on the adjacent the intervals. This requires 



9n+l . O-n+l ~ „ °n+l fa n+l 

1 c 2 e 

«n+l 

g - + ^(v n h _ 1 -V n h + ^h n )/(e-^ 
^ + ^ (V* - + /(V^^ 

«n+l e V a n+l , 







On 


a n 

— c 2 

£ 


+ t( 


Vti 




'v n h - 
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1.3. CONSTRUCTION OF A DIFFERENCE SCHEME 



Let s " = h^[ 1 - e 



Then the scheme takes the form: 



9n + ( Vjj - Vtl _ 9n\ / = 

a n \ K a n J J 

9n+l . (Vn+l- V n 9n+l\ „°n+i h n+l / 

= + 7 e * / s n+1 . 

O-n+1 \ ftn+l a-n+i I J 

Rewrite the scheme in the form: 

9n ~ Vn-1 _ 9n _ ffn+1 ( ~ Vn _ 9n+l \ g - 

s n a n fln+l \ s n+lh n +i s n+l a n+l ) 



So, 



Vj? +1 -V£ Vn-Vn-1 9n+l _^+l^±i 
1 e " 1 = e 

<?n 5n <?n+l 



S n Cl n Oi n O.J1+1 

Supplement this scheme by an approximation of the boundary conditions: 

V h = A, eV' + f(V(L)) = 0. 



In the last interval V'{x) = c^N^e^^'^ + f^, therefore, boundary 
conditions take the following form: 

V h = A, 

cga N + ^ + f(V&) = 0, i.e. 
a N 

( V N ~ V N-l) a N gNh N g N e h _ 

_a M h IL a N h N T +/l l/ ArJ 

1 - e 1 - e U7V 

So, as a result of the difference scheme has the form: 

e s = e £ — 



s n+ih n +i s n h n s n ,+ia n +i 

_ g " _|_ 3n_ _ 9n+l 

S n Cl n Oi n Q.J1+1 

V h = A, n = l,2,...,(N-l) 

"iV h JV „ a N h N „„ M V ' 

1 - e i 1 - e i ' UN 

Due to g n = giy^_i), the scheme (1.3) is nonlinear and can be linearized by 
the method of iterations [6], and solved by sweeping at each iteration. [7]. 
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1.4. EVALUATION OF CONVERGENCE 



1.4 Evaluation of convergence 

Theorem 1.4.1. Letuix) be solution of the problem (1.0.1), u h the solution 
of the scheme (1.3.2). Then there exists a constant C, independent ofe, such 
that 



T £ u = -eu" + a(x)u' + g(u), u(0) = A eu'(L) + f(u(L)) = 0, 

f e v = -ev" + a(x)v' + g(v), v(0) = A ev'(L) + f(v(L)) = 0. (1.4.1) 



We use the fact that the scheme (1.3.2) is exact for the problem (1.4.1). 
Therefore it is sufficient to estimate the proximity of these problems. 

Let z = u — v. 

Write the problem on z: 



max u n — u(x n ) < C ■ max/i, 



n n 



Proof. Rewrite the problem (1.0.1) and auxiliary problem: 



(1.4.2) 




*(0) = 0, ez'(L) + f{u{L)) - f{v{L)) = 0. 



Rewrite the problem in the form: 



-ez" + az' + g(u) - g(v) + g(u) - g(u) = 



-ez" + az' + 



g{u) - g(v) 



■ z = g(u) - g{u) 



u — V 



For z(x) we obtain the boundary problem: 



Lz = -ez" + a{x)z + g' u (9)z = g(u) - g(u), 
2(0) = 0, Dz = ez'(L) + f' u (9)z(L) = 0. 



Estimate g(u) — g(u). 
For x G A n we have: 



a(x) = a(x n -i), g(u) = g(u{x n -i)). 
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1.4. EVALUATION OF CONVERGENCE 



We obtain: 

\g(u(x)) -g(u{x))\ = \g{u(x n )) - g(u(x))\ 

\u(x n ) - u(x)\ 



< 



du 



< C \u(x n ) - u(x)\ < C 



<C \u'(s)\ ds 



u'(s) ds 



< c 



X 



•Eft — 1 



1 a(s-L) 

1 + -e~ ~ 

e 



ds < Ch n H e 2 £ 

£ 



Let /i = maXfj /t n . 
Define function ip{x): 



ifj(x) = Ch 



a(x-L) 2f3x 

e 2 £ -\- e a 



±z(x) 



Then: 

?/>(0) > 0, £>V > 0, > C/i^e^^ +C-C 2 h- Ch n ^ e ^^~ > 



By the Principle of Maximum ip(x) > 0. 
Consequently, |,s(a;)| < C/i. 



□ 



Now we have constructed a difference scheme for a nonlinear second-order 
equation with exponential boundary layer and proved its uniform conver- 
gence in the small parameter. 
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Chapter 2 



Nonlinear differential 
equations with power-law 
boundary layer 



Consider the following boundary problem: 

f (e + xf u" - f(x, u) = 0, 

l y ' \ 2.0.1) 

\u(0) = A, u(l) = B, 

under asumption that |^ > a > 0, e £ (0, 1] and the function / is continu- 
ously differentiable in its arguments. 

Let's investigate the problem (2.0.1). We will prove that the solution con- 
tains exponential boundary layer, construct a difference scheme and prove 
its convergence. 

C and Ci will everywhere designate the positive constants independent of e 
and the mesh gauge. 



2.1 Boundedness of solutions 



We prove the boundedness for the solution of the problem (2.0.1) 
Lemma 2.1.1. 

IKx)||<-||/(x 1 0)|| J 
a 

where \u{x) \ = max a |it(x)|. 
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2.2. ESTIMATION OF THE DERIVATIVE 



Proof. Write the equation (2.0.1) in the following form: 

u = f(x,0) 



T f , \2 n f(x,u)-f{u,Q) 
Lu = (e + x) u 



u 

Lu = (e + x) 2 u" — c(x) u = f{x, 0), 
where c(x) > a > 0. 

Introduce the function ip{x) = ii^fe^Ji ± u(x). 

i4(x) < 0,-0(0) > o,v(i) > o, 

Hence, according to the Principle of Maximum i/j(x) > 0. 
MiESl ± u ( x ) > o, consequently, \u(x)\ < \ ||/(x,0)||. 



□ 



2.2 Estimation of the derivative 

Lemma 2.2.1. 



\u'(x)\ < 



C 



e + x 



Proof. From the equation (2.0.1) we conclude that 

C 



\u"(x)\ < 
u"(x) = 



(e + x) 2 
f{x,u) 
(e + x) 2 ' 



Integrate (2.2.1) from £ to x: 



consequently, 



\u'{x)-u'{i)\ 
i) Let £ < x. Then 



x f(x,u(x)) 
(e + x) 2 





" -C " 


X 


< 


= C 




e + x 





e + £ e + x 



(2.2.1) 
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2.2. ESTIMATION OF THE DERIVATIVE 



ii) Let £ > x. Then 



C 



(e + x) 2 



dx, 



consequently, 



/(x)-u'(Ol <-C 



1 1 



e + £ e + x 



Thus for any £ 



Let e be so that: 



Then 



Since 



then 



\u'(x)-u'(0\ < C 



1 1 



£ + £ £ + X 



it/CO = u(i) - u (\ 



< Co and -<£<!. 



|u'(e)| <C + C 



1 1 



£ + £ e + x 



1 1 £ + X 



We took into account that 



1 1 1 

< -r < t = 2 



£ + £ " e ' 1 " 1 



2 2 



< — - < 2, 

e + £ 

|«'(0| <C + 2 + -^— < 



£ + X £ + X' 



Hence the lemma follows. 



□ 



According to Lemma 2.2.1, at the boundary x = we have the power 
boundary layer. 
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2.3. CONSTRUCTION OF THE SCHEME 



2.3 Construction of the scheme 

Let's proceed to the problem with piecewise constant coefficients. 

f(e + s)V'-/(z,«)=0, f231) 
\v(0)=A, v(l)=B, 

where f(x,u) = f n = f(x n -i,V(x ni )), at x G A n = [x n _i,a; n ]. 
To construct the scheme of equation (2.3.1) let's express v(x) as: 



v"{x) 



(e + x 



i2' 



e + x 

u(z) = -/„ • ln(e + x) + Cix + C; 



2 



Let's label v{x n -i) = V^_ x , v(x n ) = V%. 
Let's find C\ from the boundary conditions: 

j -f n - In(e + x„_i) + Cix n _i + C 2 = V^_ 1? 
\ - /„ • In(e + x n ) + C lXn + C 2 = V*, 

Cl h n -f n In ^^ = V n h -V n h _ 1 



Then, 



^ fn i £ + X n Vn Vn-l 

C\ = — In h 



fon ^ "I - X n —i h n 

We match the derivatives of solutions to the ends of adjacent intervals: 
lim V(x) = lim V'(x) 

X— >X n - X^tX n+0 

fn fn 



+ c x = — J ^±^ + c 1 



e + x n e + x n+ i 

fn fn 



Z | X lb' 



n 



x ln ( £ + Xn ) + 

^£ + X n -i) 

Vn ~ jjv-1 fn+1 fn+1 

E — T 



+i 
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2.4. JUSTIFICATION OF CONVERGENCE 



Thus we obtain: 



fn 



V. 



h 



V h i 

* 71—1 



yh 
v n+l 



V 1 



fn 



In 



£ H~ 



h 



n+l 



fn-\ 



^ fn+i ln e + x n+1 p 



yh 



e + x n _i e + x n+ \ ' h n+ i e + x n 
B, n = l,2,3,...,(iV-l); f n = f(x n . u Vti) 



2.4 Justification of convergence 



Define the mesh Vt which densities in the boundary layer in such a way that 



\e+x n -i J N " 

Let x n = A (-^f), then 



ln(e + A(- 



n 



In e + A 



n — 1 
N 



C_ 
N' 



Consequently, 



A' (f ) 1 _ C 



+ A(f) N N' 



That is 



C, A(0) = 0,A(1) = 1, 



Consequently, 



A(t) = e 



1 

1 + - 

e 



Thus, 



n 



X n — £ 



i + - 1 



n = 0, 1, 



Theorem 2.4.1. 



V h - [u} ( 



C , / 1 



Proof. Let z = u — v, where u is the solution of problem (2.0.1). Since 
the scheme (2.3.2) is exact for the problem (2.3.1) by its construction, it is 
sufficient to estimate \\v — v\\. 
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2.4. JUSTIFICATION OF CONVERGENCE 



Problem for z has the form: 

\2 II rl 



( e + x y z » _ f'(6(x))z = f(x, v) - f(x, v), 
z(0) = 0, 0(1) = 0, 



(2.4.1) 



We estimate the right side: 

f{x,v) - f(x,v) = \f(x,v) - f(x n -i,v n -i)\ < 
<C\x- Xn-tl + d \v(x) - < Ch n + Ci I 

Since \v'(x)\ < 
f{x,v) - f(x,v) 



< Ch n + Ci I — — ds < Ch n + dC 2 In ' 



^n — 1 



x + e 



e + x n _i 



In accordance with the construction of the mesh: 

f(x,v) - f(x,v) 



<Ch n + ^\n(l + - 

N \ £ 



We estimate h n . 
i) From 



Xn — £ 



we conclude that 

hjy = xn — xn-i = 1 — e 



1 + - I 



1\ * 
1 + - 1 



In(l + i) 



ii) h n increases. Therefore it is sufficient to estimate hjy. 
Then |/(s,i;)-/(x,u)| < f ln(l + ±). 
Let ^(x) = # ln(l + i) ±z(x). 

V>(0) > 0,V(1) > 0,L^(x) < 0. 

Consequently, 

|*(x)|<I.£]n(l + ± 
a N V e 



Thus, 



(7/1 
|u(x) — u(x) | < — In I 1 H — 



Since [V]q coincides with V h , we obtain the theorem. 



□ 
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Results 



At the result of the work, a first-order accurate difference scheme for a 
nonlinear equation with exponential boundary layer and nonlinear boundary 
conditions is constructed and its uniform convergence with respect to the 
small parameter is proven; a first-order accurate difference scheme for a 
nonlinear equation with power-law boundary layer and boundary conditions 
of first kind is constructed on a special mesh, and its convergence is proven 
as well. 
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